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ABSTRACT 

The presence of dust at high redshift requires efficient condensation of grains in SN 
ejecta, in accordance with current theoretical models. Yet, observations of the few 
well studied SNe and SN remnants imply condensation efficiencies which are about 
two orders of magnitude smaller. Motivated by this tension, we have (i) revisited the 
model of Todini & Ferrara (2001) for dust formation in the ejecta of core collapse SNe 
and (ii) followed, for the first time, the evolution of newly condensed grains from the 
time of formation to their survival - through the passage of the reverse shock - in the 
SN remnant. We find that 0.1 - 0.6 M Q of dust form in the ejecta of 12 - 40 M Q stellar 
progenitors. Depending on the density of the surrounding ISM, between 2-20% of 
the initial dust mass survives the passage of the reverse shock, on time-scales of about 
4 — 8 x 10 4 yr from the stellar explosion. Sputtering by the hot gas induces a shift of the 
dust size distribution towards smaller grains. The resulting dust extinction curve shows 
a good agreement with that derived by observations of a reddened QSO at z = 6.2. 
Stochastic heating of small grains leads to a wide distribution of dust temperatures. 
This supports the idea that large amounts (~ 0.1M Q ) of cold dust (T ~ 40K) can be 
present in SN remnants, without being in conflict with the observed IR emission. 

Key words: dust, extinction - shock waves - supernova remnants - supernovae: 
individual (Cassiopeia A) 



In the last few years, mm and submm observations of sam- 
ples of 5 < z < 6.4 quasars have provided a powerful way of 
probing the very existence and properties of dust within 1 
Gyr of the Big Bang. The inferred far-IR luminosities are 
consistent with thermal emission from warm dust (T < 



100 K), with dust masses ^ 10 s M Q (Bertoldi et al 
iRobson et all |2004 iBeelen et all 12006 : iHines et al.l 



2003; 



20061 ). 



Despite the uncertainties due to the poorly constrained dust 
temperatures and absorption coefficients, the estimated dust 
masses are huge, implying a high abundance of heavy ele- 
ments at z m 6, consistent with the super-solar metallicities 



these systems ( 


Pentericci et al.ll2002l 


Maiolino et al. 


2003r). 



Although high redshift quasars are extreme and rare 
objects, hardly representative of the dominant star forming 
galaxies, the above observations show that early star forma- 
tion leads to rapid enrichment of the Interstellar Medium 
(ISM) with metals and dust. 

It is difficult for the dust to have originated from low- 
mass evolved stars at z > 5 as their evolutionary timescales 
(10 s to 10 9 y r) are comparable to the age of the Univers e 
at that time (|Morgan fc Edmunds! [20031 ; lMarchenkoll2006h . 



Thus, if the observed dust at z > 5 is the product of stellar 
processes, grain condensation in supernova (SN) ejecta pro- 
vides the only viable explanation for its existence. This sce- 
nario has recently been tested through the ob servation of the 
redde ned quasar SDSSJ1048+46 at z = 6.2 (|Maiolino et all 
120041 ). The inferred extinction curve of the dust responsible 
for the reddening is different with res pect to that observe d 
at z < 4 (Small Magellanic Cloud-like. iHopkins et al.ll2004h . 
and it shows a very good agreement with the extinction 
curve predicted for dust formed in SN ejecta. 

Theoretical models, based on classical nucleation the- 
ory, predict that a few hundred days after the explo- 
sions silicate and carbon grains can form in expanding SN 
ejecta, with condensation efficiencies in th e range 0.1-0.3 
llKozasa et al l Il99ll : iTodini fc Ferrara! l200ll ; IClavton et all 
l200lT l. 

Direct observational evidences for dust production have 
been colle cted only for a lim ited numb er of SNe, such 
as 19 87A JWooden et alj Il993l'> 1999em (|Elmhamdi et all 
l2003h . and 2003ed ijSueerman et a"i1l2006l ). With the excep- 
tion of 2003gd, the dust masses derived from the IR emission 
are ~ 1O~ 3 M0, corresponding to condensation efficiencies 
which are two orders of magnitude smaller than what the- 
ory predicts. A fraction of dust could escape detection if 
it is cold and concentrated in clumps. This has been con- 
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firmed to be the case for SN 2003gd where a radiative trans- 
fer code has been used to simultaneosly fit the optical ex- 
tinction and IR emission, leading to an estimated dust mass 
of 2 x 1Q~ 2 Mq rtSugerman et al.ll2006T ). However, when ap- 
plied to SN 1987A, the same numerical model gives dust 
mass estimates whic h do not differ signifi cantly from previ- 
ous analytic results (jErcolano et al.|[2007T ). 

Similar low dust masses have been inferred from in- 
frared observations of galactic SN remnants with Spitzer 
and ISO satellite s ijHines et all 120041 : IKrause et all 120041 : 
iGreen et al.ll2004h . The consistent picture that emerges is 
that the mid- and far-IR excess observed is due to emission 
from small amounts of warm dust, with indicative tempera- 
tures T ~ 80 — 270 K and masses 3 x 10 -3 - 10~ 5 M© for Cas 
A, and temperatures T ~ 50 K and masses 3x 10~ 3 — 0.02Mq 
for the Crab nebula. Cold dust has also been detected 
through far-IR and submm observatio ns of these remnants 
jPunne et al.ll2003l ; IKrause et al.ll2004 ). However, the inter- 
pretation of these data is complicated by the strong contam- 
ination from cold dust along the line of sight, providing so 
far only upper limits of O.2M0 on the amount of cold dust 
associated to the SN remnants. 

The aim of the p resent paper is to critica lly assess the 
model developed by iTodini fc Ferraral (|200lj ) exploring a 
wider range of initial conditions and model assumptions. We 
then follow the evolution of dust condensed in SN ejecta on 
longer timescales with respect to previous theoretical mod- 
els. In particular, we are interested in understanding how the 
passage of the reverse shock affects the newly formed grain 
size distributions and masses, so as to make predictions for 
the expected dust properties from the time of formation in 
the ejecta (a few hundred days after the explosion) to its 
survival in the SN remnant, hundreds of years later. So far 
this process has received little attention, most of the stud- 
ies being dedicated to the destru ction of ISM dust grains 
caused by the SN forward shock llDraine fc Salpeterl 1 19791 : 
Ijones et al] 19941 ; Nozawa et al.ll2006T ). with the notable ex- 
ception of Dwek l|2005f ). who, on the basis of timescale con- 
siderations, finds that the reverse shock is able to destroy 
much of the initially formed dust. 

The paper is organised as follows: Sect. [2] revisits the 
dust formation models based on the nucleation theory; 
Sect. [3] describes the model adopted for the propagation of 
the reverse shock into the ejecta and shows the effect of 
sputtering on the grain size distribution and total mass; in 
Sect.|3]we compare the extinction and emission properties of 
the surviving SN dust with observations. Finally, the results 
are summarised in Sect. [5] 



2 SN DUST FORMATION REVISITED 

Models of dust formation in the ejecta of core collapse SNe 
typically predict that large masses of dust (0.1 — l.OM©) are 
formed within 1000 days from the onset of the explosion, 
when the ejecta ar e still compact (radius of order 10 16 cm; 
Kozasa et all Il99ll : ITodini fc Ferraral l200ll : iNozawa et"aH 
20031 ). If these freshly formed dust grains were distributed 
homogeneously within the ejecta, their opacity would be 
very high, with center-to-edge optical depths of order 10 2 — 
10 4 in optical wavelengths, depending on the grain material 



and size distribution. The ej ecta would thus be opaque to 
radiation produced within it (jKozasa et al.|[l99l[ ). 

Observations of recent SNe, instead, reveal extinc- 
tions smaller than a couple of ma gnitudes, which imply 
dust masses of only 10~ 4 - 10" 2 M Q l|Sugerman et al.ll2006l ; 
lErcolano et al,ll2007f ). The dust mass derived from extinc- 
tion measures could be underestimated if grains are dis- 
tributed in clumps with a small volume filling factor: for 
a given amount of grains, a clumpy distribution would pro- 
duce a lower effective extinction. However, the comparison 
between observations of dust extinction/emission and radia- 
tive transfer models shows that the neglect of clumping can 
only produce a moderate underest imation of the dust mass 
in the ejecta (|Ercolano et al.l [20071 ). 

To check whether the dust production in SNe is overesti- 
mated , we have reconsidered the model of ITodini fc Ferraral 
(2001). In the model, dust formation is investigated in the 
framework of standard nucleation theory: when a gas be- 
comes supersaturated, particles (monomers) aggregate in 
a seed cluster wh ich subsequently grows by accretion of 
other monomers l|Feder et all Il966h - For grain materials 
whose molecules are not prese nt in the gas phase, t he ke y 
species approach is adopted l|Kozasa fc Hasegawal ll987T ). 
Six materials where considered in the original work: amor- 
phous carbon (AC), iron, corundum (AI2O3), magnetite 
(Fe 3 4 ), enstatite (MgSi0 3 ) and forsterite (Mg 2 Si0 4 ). Fol- 
lowing |S^lrneid*r_et_alJ (|2004r ). we have added the formation 
of Si02 grains. SiC grains, found in meteorites and consid- 
ered to be of SN origin f rom their anomalous isotopic ratios 
l|Clavton fc Nittlerl l2004), are not considered since their for- 
mation is_jmrjed^d_by_the formation of AC and Si-bearing 
grains (INozawa et alj|2003l ). The model of lTodini fc Ferraral 
also considers the formation and destruction of SiO and CO 
molecules: while the first is necessary to study the formation 
of Si-bearing grains, the second may be a sink for carbon 
atoms that otherwise would accrete on grains. 

The ejecta are taken to have a uniform composition 
and density, with initial temperature and density chosen to 
match the observations of SN1987A. The initial composition 
depends on the metallicity and mass of the progenitor star, 
M s tar, while the dynamic is given by the mass of the ejecta 
M e je and the kinetic energy o f the explosion E^^. the models 
of IWooslev fc Weaverj jt995) were used. 

In the models of Todini fc Ferraral . the gas becomes su- 
persaturated after a few hundred days from the explosion. 
The nucleation process starts at temperature between 1800K 
(for AC) and 1200K (for Si-bearing materials). At the be- 
ginning the gas is moderately supersaturated and large seed 
clusters, made of J\f monomers, tend to form. However, their 
formation rate per unit volume (the nucleation current) is 
small. As the volume of the ejecta increases, the supersat- 
uration rate grows and smaller clusters aggregate with a 
larger formation rate. This occurs until the gas becomes 
sufficiently rarified (because of expansion and/or exhaus- 
tion of monomers in the gas phase) and the formation rate 
drops. The nucleation process, together with accretion, re- 
sults in a typical log-normal grain size distribu tions (see, 
e.g., ITodini fc Ferrarall200ll ; INozawa et al.ll2003r ). 

For materials apart from AC, the supersaturation rate 
increases quickly during the ejecta expansion, and the seed 
clusters can become very small. In ITodini fc Ferraral (|200ll ) 
seed clusters were allowed to be of any size. In this paper we 
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100 

grain size [A] 

Figure 1. Size distribution for grains formed in the ejecta of a 
SN with a progenitor star of solar metallicity and mass M s t ar = 
20M© (jV = 2). The distributions of AI2O3 and Mg 2 Si0 4 are 
shown with dashed lines for ease of identification. 



consider only clusters with M ^ 2, and introduce discrete 
accretion of monomers. While these two (more physical) re- 
quirements have a limited effect on AC grains, they alter 
the size distributions and masses of grains composed by the 
other materials. In Fig. [1] we show the size distribution of 
grains formed in the ejecta of a SN with a progenitor star of 
solar metallicity and M sta r = 20Mq . Only AC grains retain 
the usual log-normal distribution. Instead, the size distribu- 
tion of grains of other mat erials lacks the low-radiu s tail. 
Compared to the results of iTodini fc Ferraral l|200ll) . their 
total number is reduced (since larger seed clusters have a 
smaller formation rate) and their mean size is larger (since 
the monomers not allowed to form the smaller clusters are 
now available to accrete on the larger). 

It is to be noted, however, that the use of the standard 
nucleation theory is qu estionable when clusters are made 
by A/" < 10 monomers l|Draine|[l979l ; iGail et al.lll984l ). To 
check what influence this limit has on the results, we have 
run models in which the formation of clusters with M < 10 
is suppressed. The resulting size distributions confirm the 
same trend: less non-AC grains form, and of larger mean 
size. Again AC is unaffected. 

In Fig. [5] we show Md us t, the mass of dust formed in the 
ejecta of SNe of solar metallicity, as a function of M star . The 
solid line refer t o the models with J\f ^ 2. Though reduced 
with respect to ITodini fc Ferraral . still considerable masses 
of dust are formed, predominantly of AC and Fe304 grains. 
If A/star ^ 25A#0 , all the available carbon condenses in dust 
grains. In the more massive models, roughly equal amounts 
of carbon goes in grains and in CO, since the molecule 
destruction mechanism provided by 56 Co decay is reduced 
because of its low yield in the ejecta. Results are similar 
(within a factor of two) if the metallicity of the progenitor 
stars is below solar. The only distinction is the model with 
zero metallicity, where stars with M st ar 35Mq produce no 
dust (Schneider et al. 20(3). No substantial differences are 
found if a different thermal history of the ejecta is assumed: 
Mdust is still of the same order of magnitude i f densities and 
tempe ratures follow the evolution adopted bv lNozawa et al] 
(2003). As already seen in Fig. [T] imposing J\f ^ 10 results 
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Figure 2. Mass of dust formed in the ejecta of a SN as a func- 
tion of the mass of the progenitor star, for models with different 
minimum cluster size and sticking coefficient. The metallicity of 
the progenitors is solar. 



in a great reduction of the number of non-AC grains: the 
dust mass in these models is entirely due to AC, which is 
unaffected by the limit (Fig. [21 long-dashed line). 

Dust formation models depend strongly on the sticking 
coefficient a. In most of the published models, and in the 
results presented so far, it is assumed that all gas particles 
colliding on a grain will stick to it (a = 1). However, theory 
predicts that a depends on the impact energy, on the grain 
internal energy, and on the material involved: for the gas 
temper ature at which most grains form a is significantly re- 
duced (|Leitch-Devlin fc Williamslll985r i. Indeed, laboratory 
experiments on the formation of cos mic dust an alogs shows 
that a ~ 0.1 for Si-bearing grains (Gail 2003). Thus, we 
have also run models assuming a = 0.1 for all the species 
considered. By reducing a, monomers stay in the gas phase 
longer and dust formation is delayed to times when super- 
saturation is larger: typically, smaller seed clusters form. For 
M > 2, the number of non-AC grains is further reduced 
and their mass becomes negligible compared to that of AC. 
Again all available carbon is locked in AC grains, but their 
size distribution is shifted towards lower radii and seed clus- 
ters form with N < 10. For J\f > 10 (Fig. [2] dashed line) 
their mass reduces and the size distribution becomes similar 
to those of non-AC grains in Fig. [T] At least for low mass 
progenitors, the predicted Mdust are closer to the values in- 
ferred by observations. 

Clearly, the thermodynamic properties of the ejecta are 
at the limits of applicability of classical nucleation theory. 
A different approach may be needed, expecially if realistic 
a values are taken into account. Unless otherwise stated, 
in the following we will study the evolution of dust grains 
resulting from models with solar metallicity for the SN pro- 
genitors and TV ^ 2, a = 1 (so as to conform to most works 
in literature). However, we will also discuss the results for 
models with different assuptions on M and a. 
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3 SURVIVAL IN THE REVERSE SHOCK 

As the ejecta expands, a forward shock is driven into the 
ISM, which compresses and heats the ambient gas. The 
ISM becomes an hostile environment for the survival of dust 
grains preexisting the SN event, m ainly because of sputter- 
ing by collisions with gas particles |Draine fc Salpeterll 19791 : 
I Jones et alii 1994 iNozawa et al.ll2006t ). In turn, the shocked 
ambient gas drives a reverse shock in the ejecta, which, by 
about 1000 years, has swept over a considerable fraction of 
its volume. The dust within the SNe, then, has to face hos- 
tile conditions inside what had previously been its cradle. 
We study the process in this Section. 



3.1 Dynamics of the reverse shock 

iTruelove fc McKeel (Il999l ) have studied the dynamics of a 
SN remnant through its nonradiative stages, the ejecta dom- 
inated and the Sedov- Taylor. They provide analytic approxi- 
mations for velocity and position of the reverse and forward 
shocks, as a function of the kinetic energy Ekin and mass 
M e je of the ejecta, and of the ISM density prsM- We use 
here their solution for a uniform density distribution inside 
the ejecta. The values for Eki n and M e j e are the same that 
were used in the dust formation models: -Ekin = 1.2 x 10 
erg and IOM0 < M e j c < 30Mq for stellar progenitor masses 
in th e range 12-40 M B and m etallicities between zero and 
solar (|Wooslev fc Weaverlll995T ) . We study the effect of three 
different ISM environments, with pism = lO" 25 , 10~ 24 and 
10" 23 g cm -3 . 

For each model, we have divided the ejecta into N s 
spherical shells. We have assumed that all shells have the 
same initial width AR = R c j c /N s , with R c j c the initial radius 
of the ejecta. The mass of each shell is conserved throughout 
the evolution. For the j shell (counting shells outwards), the 
initial velocity of the gas at its inner boundary is given by 
homologous expansion, 



10 -Ekin 
3 Afeie' 



(1) 



where Rj is the initial radius of the inner shell boundary and 
v e j e is the velocity of the external boundary for ejecta of uni- 
form density. For practical purposes, we start our simulation 
at a time to (ideally, to — » 0), and we set R B j e = v e j e to- The 
results do not depend on the exact value of to, provided it 
is taken small enough (we use a value of order a few tens of 
years) . 

After setting the initial conditions, we study the evo- 
lution of the ejecta with time. At each time step, the re- 
verse shock goes inward through a single shell. Thus, at time 
U, the reverse shock has travelled inward through i shells, 
and lies at the inner boundary of shell j IS = N s — i — 1. 
Shells that have not been visited by the reverse shock (for 
^ j < j IS — 1) continue to follow homologous expansion, 
i.e. the inner and outer radii grow linearly with time, with 
velocity given by Eq. [T] Following the shell expansion (in- 
crease in the shell volume Vj), the shell gas density decreases 
as pj oc V7 ■ For an adiabatic expansion, the shell tempera- 



ture scales as Tj oc Vj 1 , with 7 = 5/3 jTruelove fc McKeel 
1 19991 ). Since the shock is strong, the results are independent 
on the initial choice for the gas temperature in the ejecta. 
For the shell j = j IS that has been swept over by the 
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Figure 3. Velocity, density and temperature of the ejecta at 
t = t c h as a function of the radius, for the model with progen- 
itor star of solar metallicity and mass M s tar = 20Mq (M e j e = 
I8M0) expanding in a ISM with density pism = 10~ 24 g 
cm" 3 . All quantities are nor malised to their characteristic val- 
ues ijTruelove fc McKeel Il999l l . For the model shown here = 
5800 yr, R ch = 10.7 pc, T ch = 5.7 x 10 7 K, v ch = 1800 km s _1 . 
It is also p c h = pism- The contact discontinuity marks the border 
between the shocked ejecta and the ISM swept by the forward 
shock. 



shock at time ti, we apply the standard Rankine-Hugoniot 
jump conditions for a strong adiabatic shock. The density, 
velocity and temperature change as 

7+1 , 



Ti = 2 



7 + 1 
7 — 1 r. 



'(7 + I) 2 k Urs ' 

where p'j and are the density and velocity before the shock 
(i.e. following the same evolution as for shells with j < j IS ), 
v IS is the velocity of the reverse shock in the reference frame 
of the unshocked ejecta (provided by ITruelove fc McKeel 
1999), m is the mean particle mass and k the Boltzmann's 
constant. To ensure mass conservation, the volume of shell 
j — j IS is reduced by a factor (7 — l)/(7 + 1). 

For the shells j TS < j < N IS shocked at earlier times 
t < ti, we compute the velocity Vj by interpolating between 
the velocity of the j = j TS shell and the velocity of the for- 
ward shock (in the ambient rest frame), as a function of the 
logarithm of the shell inner radius. Ve locity and position 
of the forward shock are also given by ITruelove fc McKeel 
(1999). As for shells with j < j TS , the evolution of density 
and temperature is derived from the condition of adiabatic 
expansion and conservation of the shell mass. The typical 
trends for velocity, density and temperature around the re- 
verse shock are shown in Fig. [3] 

We checked the results (in particular the assumption 
for the evolution of j > j IS shells) with the 1-D hydrody- 
namic al models of SN blast waves of Ivan der Swaluw et al.l 
l|200ll ) and with simulations kindly provided by L. Del Zanna 
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(based on the code described in iDel Zanna et al.1 [2003). 
However crude, our approximation provide a simple and fast 
solution for the density and temperature evolution of the 
ejecta during the passage of the reverse shock. Chosing an 
adequate number of shells (we use N B = 400), it agrees with 
the complete hydrodynamical solution within a factor of 2. 



3.2 Dust grain survival 

We assume that dust grains are distributed uniformly within 
the ejecta, and that the size distribution is the same every- 
where. In the shells that have been visited by the reverse 
shock, dust grains are bathed in a gas heated to high tem- 
perature (of order 10 7 — 10 8 K for the cases studied here). 
Also, the gas is slowed down and dust grains decouple from 
it, attaining a velocity relative to the gas 

2 _ 
Vdj = — r-r«rs. 
7+1 

Gas particles thus impact on dust grains transferring ther- 
mal and kinetic energy, which are of the same order of mag- 
nitude (both depending on the reverse shock velocity v TB , 
which is of order 10 3 km s -1 ). Thermal and non-thermal 
sputtering result, which erode the dust grain, reducing its 
size. Eventually, the gas drag due to direct and Coulomb 
collisions slows the grain and non-thermal sputtering weak- 
ens. In this work we consider both thermal and non-thermal 
sputtering, but we neglect the gas drag and the grain charge: 
once passed through the reverse shock, the grain retains its 
velocity relative to the gas. We can thus provide upper limits 
on the influence of non-thermal sputtering. 

The number of atoms that are sputtered off a dust 
grain per unit time is given by the sputtering rate dN/dt, 
a complex function of the gas density, temperature and of 
the nature of the dust /gas (target/projectiles) interaction 
(full expressions for dN /dt can be found elsewere, see e.g. 
iBianchi fc Ferraral lioblil ). The sputtering rate depends on 
the sputtering yield, Y , the fraction of atoms that leave the 
target per projectile collision, which is a function of the en- 
er gy of the impa c t. We use here the Y functions described 
in INozawa et al.l l|2006l ). and we consider collisions of dust 
grains with H, He and O atoms in the ejecta. The grain 
radius decreases with sputtering as 



dt 3qa 2 dt ' W 

where, q is the number of atoms in a molecule of the grain 
material, and a m is the molecule radius, computed from 
the material density and the molecule mass. The values for 
q m can be derived e asily from the acy values of Table 2 in 
INozawa et al.l {2003). At each time step, we reduce the grain 
size according to Eq.[2]in all shells that have been swept by 
the reverse shock. We follow the evolution until the reverse 
shock arrives near the center of t he ejecta: this is the l imit o f 
validity of the approximations in lTruelove fc McKed l|l999l ). 
After that, we simply assume that the ejecta expands adi- 
abatically, and we end the simulations when the sputtering 
rate becomes negligible. Since we do not include gas drag 
and grain charge, grains do not attain differential veloci- 
ties for different sizes. Thus, we have neglected destruction 
due to grain-grain collisions. However, sputtering dominates 
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Figure 4. Changes in the size distribution of AC and Fe304 
grains. For each material, the thick line is the initial size distribu- 
tion (the same as in Fig. [TJ. The thin line is the size distribution 
after the passage of the reverse shock through the ejecta. 



over this process fo r the high shock velocities considered here 
l| Jones et al.lll994h . 

Dust grains in the ionized shocked gas are heated mainly 
by collisions with electrons. If the grains are small, heating 
is stochastic and an equilibrium temperature does not exist. 
Instead, a broad temperature distribution P{T&) establishes, 
peaki ng at low te mperature but extending also to high val- 
ues (iDwekl fl98rt ). The temperature may be so high that 
dust grains sublimate l|Guhathakurta fc Draindll989l ). For 
the cases studied here, however, sublimation is negligible. 
Details on the calculation are presented in Appendix |A1 

In Fig. U] we show the initial (thick lines) and final 
(thin lines) size distributions for AC and Fe3C>4 grains in the 
ejecta of a star with M s t ar = 20M© expanding in a medium 
with pism = 10~ 24 g cm -3 . As it is evident for the (initially) 
more peaked size distributions of magnetite, sputtering pro- 
duces a leaking towards smaller sizes. The evolution of the 
size distribution is ana logous to that of ISM grains destroyed 
by the forward shock l|Nozawa et al.ll2006l ). 

In Fig. [5]we show the mass of dust that survives the pas- 
sage of the reverse shock. For the reference model (dashed 
line) , the erosion caused by sputtering reduces the dust mass 
to about 7% of its initial value, almost independently of the 
stellar progenitor model. Most of the dust (about 70% in 
mass) is consumed within one characteristic time t c h from 
the explosion, when 95% of the original volume of the ejecta 
has been reached by the r everse shock (t c h = 4— 8x 10 4 yr for 
the ejecta discussed here: iTruelove fc McKeell 19991 ). Dust in 
the inner shells is less affected by erosion, because the sput- 
tering rate is lower. A minor fraction of the dust mass, less 
than 10%, is consumed after the revers e shock bounces at 
the c enter of the ejecta (for t > 2.6t c h; ITruelove fc McKed 
Il999h . 

If the SN explodes in a denser ISM, the reverse shock 
would travel faster inside the ejecta and would encounter a 
gas at higher density. This increases the effect of sputtering. 
In Fig. [5] (long-dashed line) we see the fraction of dust mass 
that survives when pisin = 10 -23 g cm -3 : only about 2% of 
the dust mass survives. Conversely, for a lower density ISM, 
a larger fraction is left: for pisM = 10" 25 g cm -3 , it is 20% 
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Figure 5. Mass of dust that survives the passage of the reverse 
shock in the ejecta, as a function of the mass of the progenitor 
star and of the density of the surrounding ISM. The solid line 
shows the initial dust mass (same as the solid line in Fig.[2]l- 

(dotted line) . While the number of surviving grains changes 
with the ISM density, the shape of the size distributions 
remain similar in all cases, with the typical patterns shown 
in Fig. H 

No substancial change is observed in models where the 
dust was produced by progenitors of metallicity different 
from solar. Dust destruction is instead enhanced in models 
where a smaller sticking coefficient is adopted. If a = 0.1 
(Sect. [2]), only 10, 3 and 1% of the original dust mass sur- 
vives, respectively, for pism = 10~ 25 , 10~ 24 and 10 -23 g 
cm" 3 (compared to 20, 7 and 2% for a = 1.0). This is be- 
cause for smaller values of a, the dust distribution is made 
by grains of smaller radii, which are more easily destroyed. 



4 EXTINCTION AND EMISSION FROM SN 
DUST 

iMaiolino etafl l|2004) measured the reddening in the rest- 
frame UV spectrum of a z — 6.2 QSO and found it to be 
different from that of the SMC, typically used to dered- 
den the spectra of lower redshift QSOs. The measured red- 
de ning is instead compatib le with the extinction law from 
the lTodini fc Ferraral (|200lh SN dust model. We repeat here 
the same analysis using the updated dust formation mod- 
els of Sect. [2] and the final dist ributions after t he re verse 
shock passage of Sect. [3] As in Maiol ino et al. I i|2004l ). we 
der ive t he extin ction properties from the grain sizes using 
the lMid 's (|l90Sl ) theory for spherical dust grains and refrac- 
tive indexes for dust materials from the literature (references 
are provided in Table | A1[>. The proce dure is analogous to 
that adopted bv lHirashita et al (2005) when modelling the 
dust e xtinction from the SNe dust models of iNozawa et al.l 
|2003t ). 

In Fig. [|J] we show the results for dust formed in SNe 
from progenitors of solar metallicity. The grain size distri- 
butions from progenitors of different masses have been aver- 
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Figure 6. Extinction law for SN dust. The thick solid line is 
the extinction law for dust freshly formed in the ejecta. The thin 
solid line is the extinction law from dust processed by the re- 
verse shock. The curves are computed from the IMF-averaged 
size distributions of grains formed in SNe from solar metallicity 
progenitors (see text for details). The gray line and shaded area 
are the e xtinction law measur ed on a z = 6.2 QSO and its un- 
certainty ( Maiol ino et al,ll2004l) . The dashed line is the extinction 
law of the SMC JPeilll992Th 



aged over a stellar Initial Mass Function (IMF): we adopted 
the Salpete r IMF, but the result s do not depend heavily on 
this choice l|Maiolino et al-lfeOQ^ ). The thick solid line repre- 
sents the extinction law of dust as formed within the ejecta, 
without taking into account the grain processing caused by 
the reverse shock. The SN dust extinction law is still flatter 
than the SMC extinction law, but the agre ement with the 
observ ations (shaded area) is worse than in IMaiolino et ail 
(HI). This is mainly due to a change in the grain materials 
that contribute to extinction: apart from AC, present in both 
the old and new model, the rise at A < 2000 A was due to 
Mg2Si04 grains, with a minor contribution from Fe304. In 
the new model, Mg2Si04 contribution is insignificant, while 
Fe304 grains (larger than in the original model) cause the 
far UV rise. The bump at A » 2500 A is due to AC grains 
and it is typical of the optical properties derived from amor- 
phous car bon formed in an i nert athmosphere (the AGAR 
sample of IZubko et ail (|l996h ). 

During the passage of the reverse shock, Fe304 grains 
are consumed more effectively than AC grains. The resulting 
extinction law (thin solid line in Fig.0 becomes flatter, lead- 
ing to an excellent agreement with observations at A ^ 1600 
A. These results apply for ejecta expanding in a medium 
with pism = 10 -24 g cm -3 . There is no significant change in 
the extinction law if different ISM densities are considered, 
since the size distributions are similar in all cases (though 
the extinction at any given wavelength is smaller for higher 
Pism, because less grains survive). It is worth noting that 
grains with a < 20 A though as abundant as larger grains, 
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do not contribute to the extinction l aw because of their re - 
duced extinction cross section. As in iMaiolino et~aH l|2004 ). 
we find that if progenitors of metallicity lower than solar are 
considered, the difference in the resulting extinction laws are 
small and lie within 0.1 y-axis units from the lines plotted 
in Fig. [S] 

Extending calculations to the infrared, we have derived 
the dust emissivity. For all the IMF averaged size distribu- 
tions, the emissivity in the wavelength range 10 ^ A/ fim ^ 
1000 is rather featureless, and can be well described by a 
power law in wavelength of index -1.4 with ft(100/im) = 40 
cm 2 g _1 for models where all dust has been processed by 
the reverse shock. Emissivities for dust formed from progen- 
itors of a given mass are within 10% of the IMF averaged 
value, while the emissivity at the end of dust condensation, 
before any significant sputtering has occurred, is found to 
be about 20% lower. No significant dependence is found on 
the metallicity of the progenitor and on pism- In all cases, 
the emissivity is almost entirely due to the large AC grainf]. 

The amount of shock-heated dust in the ejecta can be 
derived from infrared observations of SN remnants. A no- 
table (and debated) case is that of Cas A, the remnant from 
an historical SN which shows infrared emission from the re- 
gion between the forward and reverse shocks. The identity of 
Cas A's progenitor is still highly debated. A star of 15-25 Mq 
that loses its hydrogen envel ope through winds (| Chevalier! 
120061 ') or binary interactions l|Young et al.l 120061 ) and then 
undergoes an energetic explosion can match all the avail- 
able observational constraints. In particular, the age and 
dynamics suggest a mass for the ejecta of 3 Mq , with about 
the same amount of gas reached by the reverse shock in the 
eject a and swept by the forwa rd shock in the surrounding 
ISM l|Truelove fc McKee|[i999l ). Given these uncertainties, 
and the dependence of the predicted dust masses on the 
stellar progenitor (see Fig. [SJ), we can only give a tentative 
estimate of the amount of dust predicted for Cas A by our 
model. An ejecta evolution compatible with observations can 
be obtained for a 12 Mq progenitor, provided we neglect the 
hydrogen mass. In such a model, «0.1 Mq of dust forms. 
By the age of the remnant (~325 yr), ~0.05 Mq survives in 
the region reached by the reverse shock, where it is heated 
by the hot gas. We also need to consider the contribution 
to emission from dust in the ISM reached by the forward 
shock. Typically, dust in the shocked ISM is exposed to a 
gas of similar density and temperat ure to those in the reverse 
shock l|van der Swaluw et alj |2001). For a standard value of 
the ISM gas-to-dust mass ratio, one would roughly expect a 
similar mass of emitting dust in the ISM. Thus, a model for 
Cas A remnant would have about 0.1 Mq of emitting dust. 

This mass appears to be more than an order of magni- 
tude larger than what could be derived fitting t he observed 
Spec tral Energy distribution (SED) of Cas A l|Hines et al.l 
120041 '). Using the emissivity predicted for SN dust, the flux in 
the wavelength range 10 ^ A//zm ^ 100 can be reasonably 
well reproduced with a single modified blackbody with tem- 



For the same reason, increasing the minimum cluster size Af 
and/or decreasing the sticking coefficient a does not affect the 
predicted extinction laws and emissivities, which are similar to 
those found for our reference model after the passage of the re- 
verse shock. 
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Figure 7. Synchrotr on-subtracted SED of dust emission in CasA. 
Data points are from Hincs et al. I 1120041 1. The solid line is a one- 
component modified blackbody fit to the data for A 100/im 
(T ss 100K, M d «! 4.0 X 1O _3 M ). The dashed line is a two- 
component fit with T « 110K, M d sa 3.0x 10- 3 M Q and T w 35K, 
M(j as O.IMq. The dotted line is the spectrum from stochastically 
heated dust in our model. See text for details. 



perature T = 100K, and a du st mass of 4 x 10~ 3 Mq (Fig. [7] 
solid linel. lffines et all (|2004f ) obtain for the cold, more mas- 
sive component a similar dust mass with T = 80K. However, 
the large uncertainties and the limited FIR coverage allow to 
fit, equally well, a two-component model with temperatures 
110 and 35K and masses, respectively, of 3 x 10~ 3 and 0.1 
Mq (Fig. [7] dashed line). Unfortunately Cas A lies on the 
line of sight of dense molecular clouds which do not allow 
a reliable estimate of the cold dust mass from observations 
at longer wavelengths in the FIR and sub-mm. Still, upper 
limits on the dust mass in the remnant are compatible with 
our model predictions l|Krause et al.|[2004T ). 

A broad span of temperatures is clearly needed for 
a reliable estimate of the dust mass in the remnant. In 
Fig. [7J we also show the SED of the shock heated dust in 
the CasA model (dotted line). Because of stochastic heat- 
ing (Appendix|X]), grains have temperatures mainly ranging 
from 10 to 100K. The SED cannot be easily modelled using 
a 2-component modified blackbody: the longer wavelength 
side could be described with a cold component of T ~ 60K, 
which would underestimate the dust mass by about a factor 
5; instead, a hot component at T » 150K would leave a sub- 
stantial residual in the fit at A < 10/xm. When comparing 
to the data for CasA, it appears that dust temperature in 
our models is overestimated. This could be due to an over- 
estimate of the dust stochastic heating, to a reduction of 
smaller grains with respect to the dust formation model, or 
to differences between the emission properties of true and 
modelled materials. However, the uncertainties in the ther- 
mal/dynamical history of the ejecta of CasA and the im- 
possibility of discriminating between ISM and ejecta dust 
emission in the spectrum prevent a more detailed analysis. 
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5 SUMMARY 

In the present work we have revisited the model of Todini 

6 Ferrara (2001) for dust formation in the ejecta of core 
collapse SNe and followed the evolution of newly condensed 
grains from the time of formation to their survival through 
the passage of the reverse shock. 

The main results can be summarized as follows: 

(i) The new features introduced in the dust formation 
model have only a minor impact on AC grains but signif- 
icantly affect other species (Si-bearing grains, AI2O3, and 
Fe304). For 12 - 40 Mq stellar progenitors with Z = Zq, the 
predicted Md us t ranges between 0.1 - 0.6 Mq; comparable 
values (within a factor 2) are found if the progenitors have 
Z < Zq. The dominant grain species are AC and Fe203. 

(ii) We identify the most critical parameters to be the 
minimum number of monomers, AT, which define a critical 
seed cluster, and the value of the sticking coefficient, a. As- 
suming AT ^ 10 (below which the application of standard 
nucleation theory is questionable) results in a great reduc- 
tion of non-AC grains because these species nucleate when 
the gas in the ejecta is highly super-saturated and smaller 
seed clusters form. This effect is further enhanced if a < 1: 
for a = 0.1 and stellar progenitor masses M star < 20 Mq, the 
total mass of dust is reduced to values in the range 0.001- 
0.1 Mq, comparable to those inferred from the IR emission 
at 400-700 days after the explosion for 1987A and 2003gd, 
the only two core-collapse SNe for which these data were 
available. 

(iii) Using a semi-analytical model to describe the dy- 
namics of the reverse shock, we have found that thermal 
and non-thermal sputtering produce a shift of the size dis- 
tribution function towards smaller grains; the resulting dust 
mass reduction depends on the density of the surrounding 
ISM: for pism = 10- 25 ,10 _24 ,10- 23 g cm" 3 , about 20%, 
7%, and 2% (respectively) of the initial dust mass survives. 
Most of dust consumption occurs within one characteristic 
time from the explosion, about 4 — 8 x 10 4 yr for core-collapse 
SNe. Thus, the impact of the reverse shock needs to be taken 
into account when comparing model predictions with obser- 
vations of young SN remnants. 

(iv) Averaging over a Salpeter IMF, we have derived dust 
extinction and emissivity. We find that the extinction curve 
is dominated by AC and Fes04 grains with radii larger than 
20 A. As a result, it is relatively flat in the range 1500-2500A 
and then rises in the far UV. Thus, the peculiar behaviour 
of the extinction produced by SN dust, which has been suc- 
cessfully u sed to interpret obser vations of a reddened QSO 
at 2 = 6.2 (Ma iolino et aljfeoo^l ). is preserved in the present 
model, and it is further amplified by the modifications in- 
duced by the passage of the reverse shock. 

(v) Using dust emissivity predicted by the model, we can 
reproduce the observed IR flux from the young SN rem- 
nant CasA adopting a single modified black-body of tem- 
perature T = 100 K, which imp lies a mass of war m dust 
of 4 x lO~ 3 Af , consistent with iHines et all <|2004h . How- 
ever, the limited observational coverage in the FIR allows 
to equally well reproduce the data adding a cold component 
with temperature T = 35 K and dust mass of O.IMq. Ac- 
cording to our model, such a mass of dust is what would 
be produced by a single 12 Mq star that has exploded after 
losing its hydrogen envelope, a plausible candidate for the 



highly debated CasA's progenitor. Because of the stochastic 
heating of small grains by collisions with hot gas electrons, 
dust in the shocked gas is predicted to have temperatures 
ranging from 10 to 100K. 

We conclude that our study supports the idea that 
core-collapse SNe can be major dust factories. At the same 
time, it shows that our knowledge of dust condensation and 
its survival in SN ejecta still lacks to control some critical 
parameters, which prevent reliable estimates of condensa- 
tion efficiencies, especially for the less massive progenitors. 
Within these uncertainties, the model can accomodate the 
still sparse observational probes of the presence of dust in 
SN and SN remnants. 
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APPENDIX A: STOCHASTIC HEATING FROM 
ELECTRON COLLISIONS 

We have derived t he temperature distribution P( Td) follow- 
ing the method of [Guhathaku rta fc Drainel (|l989h . to which 
we refer for a more detailed description. 

We have divided the range of possible dust tempera- 
tures into Ab bins. For the case studied here, we found suf- 
ficient to define Ab = 500 bins, logarithmically spaced in 
the range 2 < Td/K< 2000. The i-bin has temperature Td,i, 
energy Ei and energy width AEi. The energy corresponding 
to each value of Td is defined as 



E(T d ) = 




C(T)dT, 



where C is the specific heat, which is derived by fitting ex- 
perimental data. We have adopted a piecewise power-law of 
the form 

C(T) — A T B . 

The fitted values for A and B are given in Table [ATI 
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Table Al. Thermal and optical properties of SNe dust materials 



Materials 


A (erg cm 3 K v ) 


B 


Specific heat 
range in T (K) 


Refs. 


Refractive index 
Refs. 


A1 2 3 


4.44 


3 


(0, 110] 


Ditmars et al. (1982) 


Koike et al. (1995, ISAS) 




1.22 10 2 


2.29 


(110, 200] 


Chase et al. (19851 


Bcgcmann et al. (T997, compact) 




7.00 10 5 
1.45 10 7 
5.44 10 7 


0.66 
0.17 



(200, 500] 
(500, 2000] 
(2000, oo) 






Fe 3 4 


22.5 


3 


(0, 10] 


Chase et al. (19851 


Mukai (19891 




7.65 


3.47 


(10, 30] 


Shepherd et al. (19851 






3.74 10 2 
1.04 10 5 
3.95 10 7 


2.32 
1.04 



(30, 80] 
(80, 300] 
(300, oo) 


Koenitzer et al. ( 19891 




MgSiOg 


9.22 


3 


(0, 20] 


Kellev (19431 


Jaeer et al. (20031 




2.01 


3.51 


(20, 50] 


Chase et al. (19851 






7.94 10 2 
1.32 10 5 
3.47 10 7 


1.98 
0.93 



(50, 130] 
(130, 400] 
(400, oo) 


KruDka et al. (19851 




Mg 2 Si0 4 


22.7 


3 


(0, 60] 


Kellev (19431 
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The function -P(Td) is computed from the transition 
matrix between the initial and final internal energy states, 
Af,i. The probability per unit time that a dust grain in the 
energy state Ei is heated to the energy state Ef (with f > i) 
by collisions with electrons is given by 



n ira AEf x 



2{E f - Ei) 



f(E f - Ei) + if E f -Ei<E+ 



AfA = < 



2E C 



-f(Ee 



| y/1 - (E f - E t )/E cH 
1 - y/i - (Ef - E l )/E cB 







otherwise. 



In the above equation, n and m are the electron number den- 
sity and mass, respectively, and f(E) is the Maxwell distri- 
bution function defined such that nf(E)AE is the number 
of electrons of energy E per unit volume. E* is the ma ximum 
energ y that an electron can transfer to a dust grain jDwekl 
Il986l l. If Ef — Ei > Ei, it is not possible to jump from stage 
i to stage / through electron collisions; if Ef — Ei < E+ 
the jump is allowed for all electrons with energy Ef — Ei 



(which will transfer their entire energy to the grain) and for 
the electrons with energy E c g (which will transfer an energy 
i? e ffC(£'eff) = Ef — Et; the function ((E) is the fraction of 
the electron energy w hich is depos ited on a dust grain when 
E > E+, as defined in |Pweklll986l1. 

As in iGuhathakurta fc Draine] (|l989l l only cooling 
terms from level f+1 to level f are considered: 



.4 



4tt 



/,/+! 



E 



E f 



™ Qabs(a, A) B\(Ta,f+i)d\, 



where B\ is the Planck function and Q a bs(a, A) the absorp- 
tion efficiency of a gr ain o f radius a. Values for Q a bs(a, A) 
were derived using the lMiei 's l|l908l l theory for spherical dust 
grains, adopting the refractive index from the references in 
Table IAH 

In Fig. IA1I we show the temperature distributions for 
AC and Fe304 dust grains of radii a=10, 50 and 200A. 
Grains are exposed to a gas with T = 10 8 K and n = 10 
electrons cm -3 , typical conditions encountered in the ejecta 
swept by the reverse shock. Smaller grains have broader tem- 
perature distributions, with a hig h temperature t ail. 

Using the sublimation rate of iDraine fc Had l|2002l l. we 
have derived a sublimation temperature T s , defined as the 
temperature necessary to completely consume a grain in 
about 20 yr, the typical time step of our simulation. For the 
grain sizes encountered in our work, we have Td > 1200K. At 
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Figure Al. Temperature distributions P(T^) for AC and FC3O4 
dust grains in a hot gas with T = 10 8 K and n = 10 electrons 
cm" 3 . P(T(j) is shown for radii a = 10 A (broader distributions), 
50 A and 200 A (distributions spanning a smaller range of Tj). 

any time in the ejecta evolution, only a negligible fraction of 
all grains would exceed that temperature. As a result, dust 
sublimation is insignificant in our models. 



